Theory of resonant photon drag in monolayer graphene 
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Photon drag current in monolayer graphene with degenerate electron gas is studied under in- 
terband excitation near the threshold of fundamental transitions. Two main mechanisms generate 
an emergence of electron current. Non-resonant drag effect (NDE) results from direct transfer of 
in-plane photon momentum q to electron and dependence of matrix elements of transitions on q. 
Resonant drag effect (RDE) originates from q-dependent selection of transitions due to a sharp 
form of the Fermi distribution in energy. The drag current essentially depends on the polarization 
of radiation and, in general, is not parallel to q. The perpendicular current component appears if 
the in-plain electric field is tilted towards q. The RDE has no smallness connected with q and exists 
in a narrow region of photon frequency uj: \hw — 2e^| < hsq, where s is the electron velocity. 

PACS numbers: 73.50.Pz, 73.50.-h, 81.05.ue 



I. INTRODUCTION 

Though the theoretical study of two-dimensional car- 
bon has a long history [1] , [2] , 3] , [|| only after experimen- 
tal evidence of existence of graphene as a stable two- 
dimensional crystal ||, this material became very 
popular. The presence of zero gap and zero electron 
mass, combined with a rather high mobility at room tem- 
perature, makes graphene an unique material for various 
fundamental and applied problems. At present graphene 
is intensively studied both theoretically and experimen- 
tally (see e.g. reviews f9l. [TojV 

The study of graphene optics (see [H||,[l2]) is stimu- 
lated by the prediction that the absorption in monolayer 
graphene should be determined by the fundamental con- 
stant a = e 2 /he [lH, and its experimental evidence 
. The investigation of coupling between photons and 
electrons in graphene attracts now an active interest of 
the community (see e.g. 0, An observation of 

amplified stimulated terahertz emission from optically 
pumped epitaxial graphene heterostructures has been re- 
ported recently However, the photoinduced cur- 
rents, namely, photon drag and photogalvanic effects in 
graphene were beyond of interest of the researchers. In 
this paper we present the theoretical analysis of these 
effects. 

The study of light pressure on solids has rather long 
history. The simplest variant is an instantaneous trans- 
mission of photon momentum to electrons. This process 
is permitted for interband transitions or in presence of 
the "third body", for example, phonons, other electrons, 
impurities. For a free particle this process is forbidden by 
conservation laws. Small value of the photon momentum 
makes Nonresonant photon Drag Effect (NDE) extremely 
weak. 

At the same time there exists a less known variant of 
this effect, namely Resonant photon Drag Effect (RDE) 



which has no weakness of usual NDE [19| , [20| , [2l| . Res- 
onance drag occurs when some partial kinetic property 
of electron gas sharply depends on electron energy. A 
small photon momentum gives an increase of the electron 
energy, that can drastically change the relaxation time. 
This leads to a significantly different contributions to the 
electron current for electrons exited along or oppositely 
to the photon direction. In 2l| the situation was stud- 
ied for interband transitions in weakly doped GaAs when 
the electron energy approaches the energy of longitudi- 
nal optical phonon. In this case electrons exited along 
the direction of photon have larger energy than electrons 
in opposite direction. Hence, their energy can exceed the 
threshold for emission of optical phonon: they quickly 
emit phonons and stop, while the opposite electrons will 
move freely till they collide with impurity. This gives rise 
to the appearance of charge flow in the direction opposite 
to the light ray. 

Here we develop another idea for RDE based on a 
sharp Fermi distribution which forbids the transitions 
below the Fermi energy tp. This idea is illustrated in 
Fig.l. Electrons are excited from the hole cone to the 
electron cone by photons with frequency lo and wave 
vector Q. The conditions for resonant transitions are 
sk + s|k — q| = uj, Hs\k\ > ep, where k is the electron 
momentum counted from the cone point, s ~ 10 8 cm/s is 
the electron velocity, and q is a projection of the wave 
vector Q of radiation to the plane of graphene. The first 
condition determines ellipse in k plane, the second lim- 
its a part of this ellipse accessible for transitions. The 
wave vector tilts the transitions towards its direction. 
Fig. 1 shows the case when the frequency is close to 
2tp. The electrons in the figure are excited from the 
right segment of the Fermi surface contour. This results 
in electron flow rightwards. Since q -C kp the RDE 
appears when the frequency is close to 2ep, namely if 
\u) — 2ep\ < sq. Inside this window the current of RDE 



FIG. 2: (Color online) Sketch of proposed experiment (see 
text for details). 



FIG. 1: (Color online) Interband phototransitions in n- 
type graphene. Left panel: diagram of transitions in the 
momentum-energy space. The hole cone is shifted in k space 
by the photon wave vector q. The transitions are permitted 
only above the Fermi level. Right panel: projection to the 
momentum plane. Filled circle represents the Fermi sea, the 
elliptic curve corresponds to the energy conservation equation 
s|k — q| + sk = u\ only momenta outside the Fermi circle are 
permitted corresponding to the right segment of the elliptic 
curve. 



has no smallness connected with q and can be estimated 
as j ~ esr na P/(Hlj), where e and s are the electron 
charge and the velocity, ira is the opacity of graphene, 
r is the transport relaxation time and P is the light in- 
tensity. Physical meaning of this estimation is evident: 
TiraP/{fujj) is the instantaneous density of exited elec- 
trons which conserve their momentum. Being multiplied 
by the current of individual electrons es, this quantity 
gives the current density. 

Below we determine both NDE and RDE for interband 
transitions in monolayer graphene with degenerate elec- 
tron gas. Due to graphene electron- hole symmetry results 
are applicable to n- and p-type graphene. In general the 
relaxation process for electrons and holes are different 
that breaks electron-hole symmetry. For concreteness, 
we consider the n-type graphene. In this case the mean 
free time of excited electrons is much longer than that of 
holes since due to different distance from the Fermi level 
holes can easier emit phonons. Thus, the contribution of 
holes will be neglected. 

Fig. 2 illustrates a possible experiment on excitation 
of the drag current in a suspended graphene sheet placed 
in (x, y) plane. Light with frequency w, wave vector Q 
{Q = ui/c) and amplitude of electric field £ illuminates 
graphene plane. We consider transitions near the cone 
singularity. In this case the current is determined by the 
projections of the electric field and the wave vector onto 
the graphene layer. [22j These quantities are E = eE = 
(fpcos/?, £ s ) and q = (1, 0)Q sin/3, where (3 is the angle 
of incidence, £ s and £ p are amplitude components of the 



electric field £ perpendicular and parallel to the incident 
plane. We ignore small modification of field caused by 
the layer. 



II. BASIC EQUATIONS 

The current of photon drag effect can be expressed via 
the probability of transition g(k) from the hole state with 
a momentum k — q to the electron state with momentum 
k and the electron velocity v(k) = sk/k as 



4e 



dk 

4^ 



v(k)r S (k), 



(1) 



where the coefficient 4 accounts for the valley and spin 
degeneracies. The dependence on the photon momentum 
results from the momentum and energy conservation laws 
and the matrix elements for transition. For simplicity we 
put below h = 1. 

The two-band Hamiltonian near the Dirac point is 



H(k) 
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sker. 



(2) 



Here <x is the vector of the Pauli matrices. The eigenval- 
ues and eigenvectors of the Hamiltonian are e± (k) = 
±sk and *±(k) = (1, ±e ,; ^ k )/\/2, where k is the po- 
lar angle of the vector k. The different signs corre- 
spond to electrons and holes. The interaction with 
the wave is determined by the matrix elements of the 
velocity Vk-ff(k) = ser between the hole and elec- 
tron states with the momenta k — q and k: v h = 
(^_(k — q)*str^+(k)), correspondingly. 
The transition probability g(k) is 



;/(k) =^\Ev+-\ 2 S(sk + s\k-c l \-u;)9(e k -c F ), (3) 



where 9(t) is the Heaviside function. The expression for 
current Eq.([T]) can be rewritten as 
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where 



— h* — h 
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(5) 



We utilized the symmetry of the tensor ay resulting to 
inclusion of the field polarization in the combinations 
e*ej + e*ei only and independence on the degree of cir- 
cular polarization. Hence, without loss of generality one 
can consider the field as linear-polarized and e as real. 

Due to the smallness of the wave vector q, as compared 
to the electron momentum, one can expand all quanti- 
ties in powers of q. Expanding by q we can write the 
argument of the delta-function as sk + s|k — q| — cj s» 
2sk — u) — sq cos 0k (we choose the direction of axis x 
along q). At the same time, q is comparable with 2sk — u) 
and we keep ourselves from subsequent expansion of the 
delta-function. 

Expanding the tensor ay, we have 

<1 



sin k (l + - cos0 k ), 
k 



cos 2 0k 



1 ■ 2 

— sin i 

k 



> k cos k , 
9 



(6) 



2TLe(a xy ) = —2 sin 0k cos 0k — — sin 0k cos(20 k ) 
From Eq.Q we obtain for components of the current 

r min(l,a) ^ T 



-2J, 
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e 2 x [-x{l - x 2 ){\ + bx) + 26(1 - x 2 )(l 
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-x 3 (l + bx) + 4bx 2 (l-x 2 )]y, (7) 



r min(l,a) T 

j y = 2J e x e y I dx — \J\ - 

to 



| - 2x(l - x 2 )(l + bx) + 2b{3x 2 - 1)}. (8) 

Here we have introduced the following notations: 
- 2 cE 2 , , , u~2e F 



Jo 



he 87rhiu 



\e\r s, t = r\k=ki 



sq 



sq 



If t is independent on the energy of electrons then the 
integration in Eq.([7]) can be done directly. The current 
has different values inside and outside the region \u) — 
2ep\ < sq. If \u — 2ep\ < sq then we have 



U = -3 Jo VI ~ « 2 ((1 " a 2 )e 2 x + (2 + a 2 )e 2 ), (9) 

4 , 

-nM 1 



Jy 



a 2 ?/ 2 e x e 



(10) 



These values represent resonant photon drag RDE. It 
remains constant if q — > 0. The value of resonant current 
is determined by Jq. For the photon flow cE 2 /8ttHoj = 



10 19 cm 2 s 
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s, J = 1.16 • 10" 6 A/cm. This 



approximately corresponds to a power of 0.1W/cm 2 for 
photons with energy O.leV. 
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FIG. 3: Resonant photon drag current in units of Jo versus 
normalized frequency (tj — 2ep)/ sq. The solid curve shows the 
longitudinal component of current j x , the field is polarized 
along the projection of the wave vector on the plane (8 = 0) 
and j y at 9 — 7r/4. The dashed curve shows j x at 9 = n/2. 



If \u — 2ep\ > sq, then there is only NDE current. It 
is proportional to q: 
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The value of NDE is significantly smaller then the RDE 
value. 

In agreement with the simple estimates the RDE has 
always the direction opposite to the direction of light 
wave vector. Its polarization dependence is explained by 
the dependence of the directional diagram of excitation: 
most of carriers are excited perpendicular to the polar- 
ization. At the same time the Fermi sea limits the tran- 
sitions by the direction of the photon wave vector. This 
circumstances together determine lower x-component of 
current if e||q in comparison with the case e 1 q and 
also the appearance of y-component of the RDE current. 

In agreement with the system symmetry, j y exists only 
if the polarization has both e x and e y components. The 
RDE current exists in a narrow window \u> — 2ep\ < sq 
which shrinks if q — >• 0. But inside this window RDE is 
much stronger than NDE so the later can be neglected 
in this window. 

The sign of x-component of NDE depends on polariza- 
tion. This contradicts to a simple assumption according 
to which the current is mainly determined by kicks which 
photons give to electrons. The origin of this difference is 
the dependence of the directional diagram on the small 
wave vector q via the parameter a^: at some polariza- 
tions electrons prefer to be excited in opposite direction 
to q. This explains the change of sign. 

Fig. 3 demonstrates the dependence of RDE current 
components on the frequency in the window \u> — 2cf\ < 
sq where RDE exists. The current vanishes at the edges 
of the window. The component j x is larger for the polar- 
ization along the y axis. The component j y appears only 



4 



Jrde/Jq 


\ i / 2 


3\ 


4/5 


6 






/e 





FIG. 4: Polarization dependence of the RDE current at oj — 
2ef; is shown by dash curve, j y is shown by solid curve. 

for tilted polarization of the light. Fig. 4 shows the de- 
pendence of j x and j y on the angle 9 between the vector 
of polarization e and the wave vector q. 

III. DISCUSSION 

We have studied the electron contribution to the pho- 
ton drag current. In fact, in the considered system 
the hole contribution also presents. The symmetry be- 
tween holes and electrons in a neutral system means that 
these contributions double. However, the result will be 
changed if to take into account the difference between 
electrons and holes caused by their different excitation 
energy: while electrons are generated near the Fermi en- 
ergy the holes appear well below the Fermi energy. This 
leads to a strong difference between the relaxation times. 



In high-mobility samples at low temperature the mo- 
mentum relaxation time near the Fermi energy is much 
greater than far from the Fermi energy. At the same 
time, quick relaxation of excited electrons (holes) to the 
Fermi energy due to electron-electron interaction (de- 
scribed by e-e relaxation time r ee ) conserves their mo- 
menta up to the moment when excitations reaches the 
temperature layer. This results in equality of holes and 
electrons contributions to the current. And vice versa, 
electron-phonon relaxation can cancel the hole contribu- 
tion if r e _ p / l -C r ee , where T e - V h is the time of energy 
relaxation due to electron-phonon collisions. Thus, the 
obtained current should be multiplied by a factor 2 in 
the case of quick e-e relaxation and be kept unchanged 
in the opposite case. We note, that when the Fermi en- 
ergy tends to zero the system becomes symmetric. 

The RDE exists in a narrow energy range Ae hsq « 
huis/c near the Fermi energy. This means that the RDE 
is visible for temperature T < Ae. For photons with 
hjj = O.leV this gives T < 3K. 

The observation of the resonant photon drag in mono- 
layer graphene is accessible to the modern experimental 
technique that allows to investigate interesting aspects of 
coupling between photons and electrons in this material. 
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